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I. INTRODUCTION 

The solutions including scattering and/or bound states of the wave equations have been 
great interest in quantum mechanical systems [1-22]. To achieve full information about 
the system under consideration one has to investigate the bound as well as scattering state 
problem. In Ref. [3], the authors have obtained the analytical scattering state solutions 
of the ^-wave Schrodinger equation for the Eckart potential. £-wave continuum states of 
the Schrodinger equation for the modified Morse potential have been studied by Wei et.al. 
[7] where they have obtained the normalized analytical radial wave functions and derived 
a corresponding calculation of phase shifts. Chen et.al. have found the exact solutions of 
scattering states for the s-wave Schrodinger equation with the Manning-Rosen potential by 
using standard method [4]. In view of a spatially one-dimensional Woods-Saxon potential, 
the scattering solutions of the Klein-Gordon equation have been obtained in terms of hyper- 
geometric functions by Rojas et.al. and they have derived the condition for the existence 
of transmission resonances [10]. In an arbitrary dimension, Chen et.al. have presented 
the properties of scattering state solutions of the Klein-Gordon equation for a Coulomb- 
like scalar plus vector potentials [16]. In Ref. [19], low-momentum scattering in the Dirac 
equation have been studied. Villalba and Greiner [17] have investigated the transmission 
resonances and supercritical states by solving the two-component Dirac equation for the cusp 
potential. In this manner, we intend to search the transmission and reflection coefficients 
and eigenvalues of the one-dimensional Schrodinger equation for the asymmetric Hulthen 
potential (ASHp). 

The "usual" Hulthen potential [23] is one of the significant exponential potential which 
behaves like a Coulomb potential for small values of spatially coordinate. The Hulthen 
potential has many application areas in physics such as atomic physics [24, 25], nuclear and 
high energy physics [26], solid state physics [27] and chemical physics [28]. In addition, 
the Hulthen potential and its various forms are used in relativistic and non-relativistic 
regions [8, 9, 12, 13, 14]. In Ref. [8], the approximate analytical scattering state solutions 
of the Schrodinger equation with the generalized Hulthen potential for any f-state have 
been obtained. Saad [9] has studied the bound states of a spin-0 particle in .D-dimensions 
and found the normalization constant in terms of incomplete Beta function. The scattering 
solutions of the Klein-Gordon equation for the general Hulthen potential have been obtained 



and transmission resonances investigated in Ref. [13]. Guo et.al. [21] have found the 
transmission resonances for a Dirac particle in the presence of the Hulthen potential in one- 
dimension. On the other hand, solutions of the bound and scattering states of the wave 
equations for the asymmetric potentials have been recently examined [18, 22, 29]. In Ref. 
[18], the authors have investigated the low- momentum scattering of a Dirac particle in the 
presence of cusp potential. In (1 + l)-dimensions, transmission resonances in the Duffin- 
Kemmer-Petiau (DKP) equation for an asymmetric cusp potential have also been obtained 
[29]. Recently, Sogut et.al. have examined the scattering and bound state solutions of the 
DKP equation in the presence of the ASHp [22]. In the present work, we study the scattering 
and bound state solutions of the one-dimensional Schrodinger equation for the ASHp [22] 



V(x) = V 



'-*)t JL z^ + 9(*) 



Jrx 



1) 



1 — qe ax ' 1 — qe bx 

where Vq is the strength of the potential and a, b, q and q are positive parameters. 8(x) is 
the Heaviside step function and for the parameters q and q hold q < 1 and q < 1. Fig. (1) 
shows dependence of the ASHp barrier on these parameters. 

The organization of the present work is as follows. In Section 2, we search the reflection 
and transmission coefficients in terms of hypergeometric functions for the ASHp barrier by 
using the form of the wave functions for x — > ±oo. In Section 3, we obtain a condition for 
extracting energy eigenvalue for the ASHp well. This condition is a transcendental equation 
which can be solved numerically. We give some numerical values of the energy eigenvalues 
for the bound states for chosen values of the potential parameters. We summarize our results 
in Section 4. 

II. REFLECTION AND TRANSMISSION COEFFICIENTS 

The one-dimensional time-independent Schrodinger equation for a particle with mass m 
moving in a potential V(x) reads 
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2 +2m(E-V(x))H(x) = 0. (2) 

Now we look for the solution of the ASHp barrier for the region x < 0. Inserting Eq. (1) 
into Eq. (2) gives 
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Using a new variable y = qe ax in Eq. (3) one obtains the following equation 
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Taking the trial wave function 

Mv) = v*0--vYf(v), (6) 

and inserting it into Eq. (4) we have 

2/(1 - y)f"(y) + [1 + 2/i - (2/x + 2v + l)y] /'(y) - (/i + zv + 7 )(/i + z/ - 7)7(2/) = , (7) 

which has the form of the hypergeometric-type equation [30] 

s(l - s)x" + [Ca " (Ci + C2 + l)*]x' - C1C2X = . (8) 

whose solution is given as 2 -^i = (Ci? C2; C3; s \ So, comparing Eq. (7) with Eq. (8) gives us 
the solution 

f(y) = A 12 F 1 (^ + z/-7,/i + zy + 7;l + 2fi;y) 
+A 2 y' 2 " 2 Fi (-// + 1/ - 7, -// + 1/ + 7; 1 - 2^; y) , (9) 

and the whole solution for the region x < 

^L(y) = Ay^i -yy 2 F 1 ( t i + v- 1 ,ix + v + rA + ^ y) 

+A 2 y-»(1 - y ) u 2 F 1 (- f i + v - 7, -/1 + j/ + 7; 1 - 2 M ; y) , (10) 

where 



/i = i-; fc = v / 2m£; 1/ = 1 ; 7 = -W2m f £+ — J . (11) 

We have to obtain the asymptotic form of the above wave function since we search the 
reflection and transmission coefficients. Asi-} —00, y — > and (1 — y) u — > 1, we obtain 
from Eq. (10) 

ij) L {x -► -00) ~ A^e "* + Aag-^e -0 ^ ~ Arf^e** + A 2 g- ifc/a e-' ifca; , (12) 



where we have used 2-^1 = (Ci, C2; C3J 0) = 1. 

To obtain the solution of the ASHp barrier for the region x > we insert Eq. (1) into 
Eq. (2) and get 
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Defining the new variable z = qe gives us 
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where 
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By using a trial wave function iPr(z) = 2 Ml (l — z) Ul h(z) in Eq. (14) we obtain the whole 
solution of the ASHp for the region x > 

iP R (z) = A 3 z^(l - z) Ul 2 F 1 (fi 1 + v 1 - 7!,^ + v x + 71; 1 + 2fn; z) 
+A 4 z-" 1 (1 - z) ui 2 F 1 (-fi 1 + 1/1 - 71, -fii + vi + 7i ; 1 - 2/ii; *) , (16) 

where 



/i 1= 4; k = V2^E; v x = 1 ; 71 = ^/2m f £ + 5 J . (17) 

In order to define a plane wave travelling from left to right we have to set A3 = in Eq. 
(16), so 

il> R (z) = A 4 z-^(1 - zY 1 2 F 1 (-fi 1 + ui - 71, -iii + fi + 7i! 1 - 2jUi; «) . (18) 

Now we give the form of the wave function at x — > +00 for region x > 0. As x — > +00, 
z — > and (1 — z) Ul — > 1, we have from Eq. (18) 



^(x ->■ +00) ~ Ai(g)-" 1 e 6 " ix ~ Ai^-^V** . 



(19) 



As a result we can summarize the wave function for the limit x — ¥ ±00 from Eq. (12) and 
Eq. (19) as 
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The wave function in Eq. (10) can be written as i/jl = i^inc + VVe/ in the limit x — > — oo 
where if)i nc is the incident and ift re f is the reflected wave. Similarly, as x — > +00 the wave 
function in Eq. (18) is ipn — 4>trans where iptrans is the transmitted wave. These definitions 
give us the reflection and transmission coefficients as 
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In order to give the explicit expressions for the coefficients used in the above equations we 
need to use the continuity conditions on the wave function given as ipn(x — 0) = iPl(% — 0) 
and ip'n(x = 0) = if)' L (x = 0) where prime denotes derivative with respect to x . The 
matching of the wave functions at x = gives 



A1C1F1 + A 2 C 2 F 2 = A A C 3 F 3 



(22) 



and the matching of derivatives of the wave functions reads 



aqA 1 C 1 (D 1 F 1 + D 4 F 4 ) + aqA 2 C 2 {D 2 F 2 + D 5 F 5 ) 

= bqA 4 C 3 {D 3 F 3 - D 6 F 6 ) 



(23) 



where we have used the property of the hypergeometric functions as 4- 2-^1 (Ci) C2; C3! s ) — 

^ 2 F 1 (Ci + i,C2 + i;Ci + i;s). 

Combining last two equations we obtain the followings for the coefficients written in Eq. 
(21) 

A 2 C\ [bqF 1 (D 3 F 3 - D G F 6 ) - aqF 3 (D 1 F 1 + D 4 F A )} 



A x C 2 [aqF 3 (D 2 F 2 + D 5 F 5 ) - bqF 2 (D 3 F 3 - D 6 F 6 )} ' 
A 4 ^ gqd [F 1 (D 2 F 2 + D 5 F 5 ) - F 2 (D 1 F 1 + D^)} 
A 1 C 3 [aqF 3 (D 2 F 2 + D b F b ) - bqF 2 (D 3 F 3 - D 6 F 6 )] " 
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where the following abbreviations in the above equations have been used 

C\ = q"(l-qr-, C 2 = «-»(! - qf ; C 3 = (s)-«(l - «T , (26) 
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D _^ (-/^i + ^1 - 7i)(-/^i + ^1 + 7i ) (27) 

6 1-2/i! ' l ' 

F 2 (n, v,7,q) = 2 F 1 (-fi + iy--f, -/i + i/ + 7; 1 - 2/x; g) , 

F 3 ( / ui, 1/1,71, g) = 2 ^i(-/Ui + 1/1-71, -fjii + 1/1 + 71; 1 - 2^i; g) , 

F 4 (ju, 1/, 7, g) = 2 Fi(n + i/ - 7 + 1, /U + i/ + 7 + 1; 2 + 2/x; g) , 

^5 (A*, z/, 7, g) = 2^i(-/i + 1/ - 7 + 1, -A* + " + 7 + 1; 2 - 2/i; g) , 

^e(A*i, ^1, 7i, 9) = 2^i(-A i i + v \ ~ 7i + 1, -A*i + ^1 + 7i + 1; 2 - 2jUi; g) . (28) 

Inserting Eqs. (26), (27) and (28) into Eq. (21) gives the explicit expressions of the trans- 
mission and reflection coefficients. Figs. (2)- (4) show different variations of these coefficients 
according to the energy and also potential strength for various potential parameter values. 
It is seen in Fig. (2) that the unitarity condition, R + T = 1, is certainly satisfied. Fig. 

(3) and left panel of Fig. (4) show that the dependence of the transmission coefficient on 
the potential parameters is very similar which means that it goes to zero for relatively lower 
values of energy while goes to unity for higher values of energy. The right panel of Fig. 

(4) displays the dependence of the transmission coefficient on the strength of the potential. 
According to this plot, the transmission probability of the particle from the barrier is exactly 
one if the height of the potential is zero as expected. This probability goes to zero with 
increasing value of the strength of the potential. 

III. ENERGY EIGENVALUES 

In this section, we deal with the bound state solutions of the ASHp well which means 
that Vo — ¥ — Vq. The equation (3) for x < turns into 
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Using the transformation y = qe ax and taking the trial wavefunction ip(y) = y' I2 (l — yY 2 9{y), 
the solution of Eq. (29) becomes 

g(y) = A 5 2 F 1 (fi 2 + v 2 ~ I2, ^2 + ^2 + 72; 1 + 2^ 2 ; v) 
+A Q y- 2 ^ 2 Fi(-H2 + "2- 72, -1x2 + z/ 2 + 72; 1 - 2^ 2 ; y) , (30) 

with the parameters 



A*2 = -V2mE = fi; v 2 = v = \\ 72 = -\2m ( £ I . (31) 

a a y V g / 

and the complete solution of Eq. (29) is given 

Mv) = A>2T(1 - vY 2 2^1(^2 + v 2 - 7 2 , // 2 + z/ 2 + 72; 1 + 2/i 2 ; y) 
+A 62 /^ 2 (1 - y)" 2 2 Fi(-M2 + ^2 - 72, -H2 + "2 + 72; 1 - 2^2; y) • (32) 

Next, we search the solutions of the following form of Eq. (13) for x > 
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By using the variable z = qe bx and putting ip(z) = z M3 (l — z) U3 u(z) in Eq. (33) we obtain 

u(z) = A 7 2 F 1 (ij,3 + v 3 - 73, fi 3 + 1/3 + 73; 1 + 2^3; z) 
+A%z~ 2 ^ 2 ^i(-At3 + ^3 - 73, -A* 3 + ^3 + 73; 1 - 2^ 3 ; «) , (34) 

with the parameters 



/i 3 = -\JlmE = //1 ; z/ 3 = z/ x = 1 ; 73 = -W2m ( £" - -^ j . (35) 

Finally, we obtain the complete bound state solution of the Schrodinger equation for x > 

ij) R (z) = A 7 z^(l - z)" 3 2 F 1 (fi 3 + iy 3 - 73, /is + ^3 + 73; 1 + 2/i 3 ; «) 
+A 8 ^ 3 (1 - z) V3 2 Fi(-/i 3 + ^3 - 73, -/is + ^3 + 73; 1 - 2/i 3 ; 2) • (36) 

In order to represent the wavefunctions in Eq. (32) and Eq. (36) of the bound state solutions 
they satisfy the boundary condition being zero at infinity which gives A e = A s = and we 
obtain 

Mv) ~ W 2 (! - VY 2 2F 1 {fi 2 + v 2 - 72, ix 2 + v 2 + 72; 1 + 2^ 2 ; y) , (37) 

i> R (z) ~ A 7 ^ 3 (l - zY 3 2^1(^3 + ^3 - 73, ^3 + ^3 + 7s; 1 + 2^ 3 ; z) . (38) 



Matching last two expressions in x = requiring continuity of the wavefunction and of 
its first derivative gives 

A 5 F 1 (fi 2 ,^2,'l2,q) - A 7 F 2 (fi 3 ,V3,'l3,q) = 0, (39a) 

4 ( — -T^—) F i(^2, »2, 72, q) + F 3 (^ 2 , ^2,72, q) 
-A 7 U^--^-\ F 2 (^ 3 ^3,73,g) + i ? 4^3^3,73,g)}=0. (39b) 

where 

Fi(fi 2 ,v2,j2,q) = q^ 2 (l -qY 2 2Fi(n2 + "2- 12,^2 + ^2 + 12; 1 + 2^2; q) , 

F 2 (V3, "3, 73, q) = (q) M3 (l - qT 3 2^1(^3 + ^3 - 73, /x 3 + z/ 3 + 73; 1 + 2^ 3 ; g) , 

r- / X ho/1 \^2 (^2 + ^2 -72) ((^2 + ^2 +72) 

F 3 {fi 2 ^2,l2,q) = q^{l-q) U2 — ^y^ 

x 2 Fi (jLt 2 + ^2 - 72 + 1, At2 + ^2 + 72 + 1; 2 + 2^ 2 ; g) , 

^4(^3,^3,73,9) = (g) M (i-g) 3 7-7-y 

x 2 -Pi(A*3 + ^-73 + l,A*3 + ^3 + 73 + l;2 + 2/i 3 ;g). (40) 

Equation (39) has a nontrivial solution only if its determinant is zero. Using this equation, 
one can determine the energy eigenvalues of the ASHp well numerically. Here, we give our 
numerical results for the energy eigenvalues as a list for some values of the parameters, for 
example, m — l,a — 0.5, b = 0.75, Vq — 5,q — 0.1 and q = 0.5 taking into account that 
- | Vb| < E < 0: E x = -2.453010, E 2 = -2.251290, E 3 = -0.924802, E A = -0.491271, E 5 = 
—0.001356 (in atomic unit). 

IV. RESULTS AND CONCLUSIONS 

We solve the one-dimensional Schrodinger equation for the asymmetric Hulthen poten- 
tial. We find the transmission and reflection coefficients for the ASHp barrier in terms of 
hypergeometric functions and give some plots showing the dependence of these coefficients 
on the potential parameters a, b, q, q, Vq and on the energy E. We observe that the unitar- 
ity condition is exactly satisfied in all cases. We also compute the energy eigenvalues for 
the bound states extracting an eigenvalue equation which can be solved numerically. We 
calculate five different energy eigenvalues by taking into account that — |Vq| < E < 0. 
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FIG. 1: plots of ASHp for different values of the potential parameters for 
Vq = l,q = q = 0.5 (left panel) and Vq — 1, a — b — 0.5 (right panel). 




FIG. 2: transmission (T) and reflection (R) coefficients varying with E for 
a = 0.4, b = 0.5, q = 0.6, q = 0.7, m = l,V = 2. 
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FIG. 3: variation of the transmission coefficient with different potential parameters for 

m = l,V = 2 (left paneha = b = 0.5; a = 0.8, b = 0.3; a = 0.3, b = 0.8; q = q = 0.7 right 

panel:? = q = 0.5; g = 0.6, g = 0.4; q = 0.4, g = 0.6; a = b = 0.5). 
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FIG. 4: variation of the transmission coefficient with energy E and potential parameter V$ 

toi a = b = q = q = 0.5, m — 1. 
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